Exact statistical properties of the Burgers equation

The one dimensional Burgers equation in the inviscid limit with white noise initial condition is revisited. The one- and two-point distributions of the Burgers field as well as the related distributions of shocks are obtained in closed analytical forms. In particular, the large distance behavior of spatial correlations of the field is determined. Since higher order distributions factorize in terms of the one and two points functions, our analysis provides an explicit and complete statistical description of this problem.Comment: 21 pages, 6 figures include

Ballistic aggregation: a solvable model of irreversible many particles dynamics

The adhesive dynamics of a one-dimensional aggregating gas of point particles is rigorously described. The infinite hierarchy of kinetic equations for the distributions of clusters of nearest neighbours is shown to be equivalent to a system of two coupled equations for a large class of initial conditions. The solution to these nonlinear equations is found by a direct construction of the relevant probability distributions in the limit of a continuous initial mass distribution. We show that those limiting distributions are identical to those of the statistics of shocks in the Burgers turbulence. The analysis relies on a mapping on a Brownian motion problem with parabolic constraints.Comment: 23 pages, 6 figures include

On the adiabatic properties of a stochastic adiabatic wall: Evolution, stationary non-equilibrium, and equilibrium states

The time evolution of the adiabatic piston problem and the consequences of its stochastic motion are investigated. The model is a one dimensional piston of mass $M$ separating two ideal fluids made of point particles with mass $m\ll M$. For infinite systems it is shown that the piston evolves very rapidly toward a stationary nonequilibrium state with non zero average velocity even if the pressures are equal but the temperatures different on both sides of the piston. For finite system it is shown that the evolution takes place in two stages: first the system evolves rather rapidly and adiabatically toward a metastable state where the pressures are equal but the temperatures different; then the evolution proceeds extremely slowly toward the equilibrium state where both the pressures and the temperatures are equal. Numerical simulations of the model are presented. The results of the microscopical approach, the thermodynamical equations and the simulations are shown to be qualitatively in good agreement.Comment: 28 pages, 10 figures include

Ballistic Annihilation

Ballistic annihilation with continuous initial velocity distributions is investigated in the framework of Boltzmann equation. The particle density and the rms velocity decay as $c=t^{-\alpha}$ and $=t^{-\beta}$, with the exponents depending on the initial velocity distribution and the spatial dimension. For instance, in one dimension for the uniform initial velocity distribution we find $\beta=0.230472...$. We also solve the Boltzmann equation for Maxwell particles and very hard particles in arbitrary spatial dimension. These solvable cases provide bounds for the decay exponents of the hard sphere gas.Comment: 4 RevTeX pages and 1 Eps figure; submitted to Phys. Rev. Let

Phase transition in a spatial Lotka-Volterra model

Spatial evolution is investigated in a simulated system of nine competing and mutating bacterium strains, which mimics the biochemical war among bacteria capable of producing two different bacteriocins (toxins) at most. Random sequential dynamics on a square lattice is governed by very symmetrical transition rules for neighborhood invasion of sensitive strains by killers, killers by resistants, and resistants by by sensitives. The community of the nine possible toxicity/resistance types undergoes a critical phase transition as the uniform transmutation rates between the types decreases below a critical value $P_c$ above which all the nine types of strain coexist with equal frequencies. Passing the critical mutation rate from above, the system collapses into one of the three topologically identical states, each consisting of three strain types. Of the three final states each accrues with equal probability and all three maintain themselves in a self-organizing polydomain structure via cyclic invasions. Our Monte Carlo simulations support that this symmetry breaking transition belongs to the universality class of the three-state Potts model.Comment: 4 page

Defensive alliances in spatial models of cyclical population interactions

As a generalization of the 3-strategy Rock-Scissors-Paper game dynamics in space, cyclical interaction models of six mutating species are studied on a square lattice, in which each species is supposed to have two dominant, two subordinated and a neutral interacting partner. Depending on their interaction topologies, these systems can be classified into four (isomorphic) groups exhibiting significantly different behaviors as a function of mutation rate. On three out of four cases three (or four) species form defensive alliances which maintain themselves in a self-organizing polydomain structure via cyclic invasions. Varying the mutation rate this mechanism results in an ordering phenomenon analogous to that of magnetic Ising model.Comment: 4 pages, 3 figure

Statistics of Largest Loops in a Random Walk

We report further findings on the size distribution of the largest neutral segments in a sequence of N randomly charged monomers [D. Ertas and Y. Kantor, Phys. Rev. E53, 846 (1996); cond-mat/9507005]. Upon mapping to one--dimensional random walks (RWs), this corresponds to finding the probability distribution for the size L of the largest segment that returns to its starting position in an N--step RW. We primarily focus on the large N, \ell = L/N << 1 limit, which exhibits an essential singularity. We establish analytical upper and lower bounds on the probability distribution, and numerically probe the distribution down to \ell \approx 0.04 (corresponding to probabilities as low as 10^{-15}) using a recursive Monte Carlo algorithm. We also investigate the possibility of singularities at \ell=1/k for integer k.Comment: 5 pages and 4 eps figures, requires RevTeX, epsf and multicol. Postscript file also available at http://cmtw.harvard.edu/~deniz/publications.htm

Aging and its Distribution in Coarsening Processes

We investigate the age distribution function P(tau,t) in prototypical one-dimensional coarsening processes. Here P(tau,t) is the probability density that in a time interval (0,t) a given site was last crossed by an interface in the coarsening process at time tau. We determine P(tau,t) analytically for two cases, the (deterministic) two-velocity ballistic annihilation process, and the (stochastic) infinite-state Potts model with zero temperature Glauber dynamics. Surprisingly, we find that in the scaling limit, P(tau,t) is identical for these two models. We also show that the average age, i. e., the average time since a site was last visited by an interface, grows linearly with the observation time t. This latter property is also found in the one-dimensional Ising model with zero temperature Glauber dynamics. We also discuss briefly the age distribution in dimension d greater than or equal to 2.Comment: 7 pages, RevTeX, 4 ps files included, to be submitted to Phys. Rev.